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We present a simple example of quantum control in Bose-Einstein condensates via Feshbach 
resonance. By tuning an initially positive scattering length to zero, it is possible to generate os- 
- - - dilatory motion of the condensate that results from quantum interference. The density oscillation 

^ID , is accompanied by a periodic enhancement of the quantum mechanical squeezing of the amplitude 

' quadrature. 
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With the rapid technological advances in quantum physics and nanoscience, it is becoming increasingly necessary 
to develop methods of quantum control based on the laws of quantum mechanics. Owing to the microscopic scale 
involved, one of the main avenues so far for the manipulation and control has been the interaction of the quantum 
I— I, system with an electromagnetic wave. Prominent examples include atom optics and cavity QED [1]. 

■ The Bose-Einstein condensate (BEC), created experimentally almost a decade ago, is a unique mesoscopic quantum 
O \ system which is being studied intensely by a great number of researchers from various disciplines. Many experimental 
^ ■ results have been reported in quick succession since the initial production of BEC [2] . One of the important, early BEC 
I experiments was that of inducing Feshbach resonances to tune the s-wave scattering length [3-9] . Recent experiments 
. have provided the opportunity to tune the scattering length from the repulsive (a > 0) through the ideal gas (a ~ 0) 
^ ' to the attractive (a < 0) regime. Feshbach resonance has become an indispcnsible tool in the studies of BEC ranging 
, from the production of molecular BEC [10,11] to the observation of a "bosenova", a phenomenon closely related to 
^ ' the supernova [12-14]. A large number of theoretical studies on the collapse and the subsequent "explosion" of BEC 
O \ has since also appeared in the literature [15]. 

I ^ I ■ A Feshbach resonance occurs when the energy of a pair of atoms in the condensate matches that of a metastable 
molecular state. The s-wave scattering length, which is a function of the energy difference between the atomic and 
the molecular states, may be varied near a Feshbach resonance using magnetic field: 
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where, for the case of ^''^Rb in the f = 2, mj — —2 hyperfine state, a Feshbach resonance characterized by the width 
^f) AB Ri 11. 6G and the off-resonant background scattering length ~ — 450ao where ao is the Bohr radius, is found at 
: Bq « 154. 9G [14]. We show in this paper that, in theory, tuning the scattering length of BEC to zero via Feshbach 

resonace (e.g. by setting the magnetic field to B = AB + Bq) constitutes a simple method of quantum control - 
j observable changes are induced in the density of the BEC, and in addition, quantum mechanical squeezing [16] is 
■ realized. 

I A recent, related theoretical work has considered periodic sign-changing modulation of the scattering length similar 

to dispersion management schemes in fiber optics [17]. This method was found to induce "breathers" which oscillate 
Q ' between the Thomas-Fermi (TF) and Gaussian configurations. Our proposed scheme, although similar, does not 
^ require changes in the sign of the scattering length. We propose a situation in which an initially positive scattering 
ILJ length a > is switched to e where e is negligibly small but positive. Because a is at all times non-negative, this 
precludes the possibility of instability and collapse introduced when a < [13,15]. We study the subsequent evolution 
of the condensate as a is kept constant at e w 0, away from the resonance. 



X 

H _ As in Ref. [17], we assume that atom losses are minimized by careful experimental control [8,10,14], and that we 
. . . 1 are in an off-resonant regime in which any molecules formed are rapidly destroyed, overdamping any atom-molecule 
population oscillations. In such a regime, the formation of a molecular condensate is suppressed and an accurate 
description of the condensate may be assumed to be provided by the Gross-Pitaevskii Equation (GPE). We write the 
time-independent GPE as follows 



2M- +^(^)+9\Mr)V 



Mr) = A'^c(r), (2) 



where fi is the chemical potential, V{r) is the harmonic confining potential of the form ^Mujfr^, and g ~ AirNah? /M 
is the interparticle coupling constant with the s-wave scattering length a and the total number of atoms N . ipc 
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represents the condensate wave function which is normalized to one, J |'^c(r)|^rfr = 1, as we have included the total 
number of condensate atoms in the coupling constant. 

If the scattering length a is set equal to zero suddenly (within a time scale t ^ 1 /cut , to avoid any adiabatic following) , 
the eigenstates of the system become the harmonic oscillator states, <^n(r), with the corresponding eigenvalues nhcot- 
It is straightforward to see that the equation obeyed by Vc(r) after the switch to a ~ can be written as a set of 
equations 

i^-^ = nTwJtan, (3) 

where an = j (f)'^{r)ipc{j^)dr are the expansion coefficients of the BEC ground state expanded using a complete set of 

eigenfunctions {0„(r)}, n = 0, 1, The coefficients a„ satisfy the normalization relation J2n ["^"P = 

arbitrarily be chosen to be real. The solution to Eq. (3), which assumes that the switch to a occurs at time t = 0, 
is trivially given by = a„(0) exp(— mcjff). This causes the condensate wave function to become time-dependent: 
ipc{^, t) = q;„(0) exp(— mwtt)0„(r). In particular, the density of the BEC ground state is now given as 

|Vc(r, t)\^ = J2 «™(0)"n(0) exp[-i(n - m)wttWm{r)Mr) (4) 

= ^|a„(0)n,^„(r)|2+ 2a„(O)a„(O)cos[(n-m)wti]0„(r)0„(r), (5) 

n n>m 

where we have used the fact that a„(0) and ^^(r) are real. Equation (5) clearly demonstates a time-dependent density 

oscillation resulting from quantum interference. The form of the density is suggestive of a diffraction grating in time, 
in which various components of the condensate wave function interfere with one another. The "visibility" V, of the 
temporal interference at the center of the trap may be calculated by defining the "intensity" / = |V'c(0,t)p: 

+ Imin E„|an(0)|2|</)„(0)|2 • 

The interference of various contributing components implies that the initial profile of the condensate is an important 

factor in determining the character of the density oscillation and the associated visibility. For a small BEC far from 
the TP limit, ipc may well be approximated by a Gaussian. This implies that amongst all the expansion coefhcients, 
ao(0) and a2(0) should be dominant. The second lowest coefficient, ai{0), that corresponds to the first excited state 
(?!)i(r) is negligible since the condensate in a harmonic trap is an even function. Assuming all the other coefficients 
are very small, we get |V'(r,<)P « |ao(0)|2|,/)o(r)|2 + |a2(O)|2|02(r)P + 2ao(0)a2(0) cos(2wtt)(/)o(r)02(r). Thus the 
condensate follows an oscillation at frequency 2ujt with the corresponding visibility at the center of the trap given 
as V = 2ao(0)a2(0)(/>o(0)(/>2(0)/[ |ao(0)|^|(/)n(0)p + |a2(0)|2|</)2(0)|2]. Figure 1 shows the maximum density variation 
of the BEC after tuning the scattering length to zero, with the top figure giving the result for a small condensate 
and the bottom figure for a condensate in the TF approximation. As expected, different oscillation patterns are 
observed whereas the frequency of oscillation is 2(jjt in both cases. This implies that in the TF approximation, the 
lower coefficients are still dominant, although with different weightings from that for a small condensate. 

It is found that the Feshbach resonance may also be used to control the amplitude and phase quadratures of 
the condensate mean field. In particular, we consider the changes induced in their respective variances i.e. the 
amount of quantum mechanical squeezing [16]. Defining the usual harmonic oscillator annihilation (creation) op- 
erator, b (b^), such that 60o(r) — 0, b4>n{''^) = V^^n-i{^) and 6Vn(r) = ■\AiH-T?^>n+i(r), we note that, using 
|^/)£,(r, t)) = a„(0) exp(— ina;ti)|(/)„(r)), the number of atoms is conserved as to be expected: {ipd^, t)\b^b\ipc{^, t)) = 
(V'c(r,0)|6t6|V'c(r,0)). On the other hand, 

(Vc(r, i)|6|Vc(r, t)) = (Vc(r, 0)|6|Vc(r, 0))e-^-**, (7) 

implying that the quadrature amplitude of the field (Xi) = ^ {b+b^) undergoes a time-dependent oscillation of the form 
{Xi{t)) = (Xi(0)) cos{u>tt) while the quadrature phase of the field is given by {X2{t)) = -^{b— b^) = {X2{0)) sm[iJtt). 
The variance (AX)^ = {X'^) — (X)"^ of these quadratures are given by 

[AXi(t)]2 = i [C{b\b) + C{b, b^) + C{b, 6)e-^2'-** + C{b\b^)e'^^^'] , (8) 



C{b^,b) + C{b,b) cos{2L0tt) + ^ 



(9) 
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[AX2(t)]2 = ^ [C{b^,b) + C{b, 6+) - C{b, 6)e-'2'^'* - 6^)6'^'"**] , 



C(&^6)-C(6, b) cos(2ojtt) + 



(10) 
(11) 



where we have introduced a shortened notation for the correlation between two operators C{X, Y) = (XY) — {X){Y), 
and used the fact that C{b, b) = C{b\b'^) for a„ reaL In terms of the normahzed expansion coefficients, 



C{b, b) = Y^ V(" + 2)(n+ l)Q!„a„+2 - ^ ^/ {n + l){m + l)a 

n n.m 

C{b\b) = ^ Vnal - ^/ {n + 1)(to + l)a„+ia„a„a„+i, 



(12) 
(13) 



and since a„ « when n is odd for the case of a harmonically trapped condensate (an even function), the correlations 

may also be written as C{b,b) » \/ (2n + 2)(2n + l)a2n«2n+2 and C(6^,6) w X^n 'V^2'^Q:2ri- Equations (8-11) 
indicate a time-dependent oscillation of the variances and hence quantum mechanical squeezing of these quadratures. 
The change from the original variances when the scattering length was not tuned is found to be 



[AX,{t)r-[AX,{0)f 



C{b,b) 



[cos{2LUtt) - 1] , 



(14) 



[AX2(i)]'-[AX2(0)]2 = 



C{b,b) 



[1 - cos{2uitt)] . 



(15) 



The maximum changes in the variances imply that the amplitude quadrature Xi is squeezed periodically at the expense 
of the phase quadrature as a result of the change in the scattering length. The variance of the phase quadrature is 
never reduced from that of the original value. Plots of these differences are presented in Fig. 2. 
The induced changes in the quantum state of BEC may be visualized using the Q-function [16,18] 



Q{0,f3*) = {P\m = e 



E 



n'. 



(16) 



where |/3) = e-l^^l'/^ 

denotes the usual coherent state with = and p = |-0c)('0c|. Figure 

3 shows the time development of the Wigner contour for the small (Gaussian) and large (TF) condensates, at times 
t = and t = it/{2uji). It is seen that the formerly stationary Wigner contour is now time-dependent and undergoes a 
rotation about its center at frequency w*, confirming the oscillatory behavior of the quandrature squeezing. It should 
also be noted that, although the density variation of Fig. 1 may be observed relatively straightforwardly, the more 
experimentally involved quantum tomography [19] or a projection synthesis technique [20] are required in order to 
completely measure the Q-function directly. 

We note that similar idea may be applied to a system of rotating BEC. In such a case, the corresponding time- 
independent GPE is given by 



2M ' 
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■5|V'rot(r)|^ 



V'rof(r) = MV'rot(r), 



(17) 



where ojr is the frequency of harmonic potential in the xy plane; r = [x, y), = \r x p| is the angular momentum, 
and f2 is the rotational frequency of the trap. The eigenstates of the system after tuning the scattering length to zero 
are then the Landau states (^^„(r): 



(r) 



g|z|V(2<)^m^ng- 



(18) 



with eigenvalues e„,m = fi-[(wr + ^)n+ (w^ — fl)m + ijJr], where z = {x + iy)/ar, = ^/fijMujy, d± = {ar/2){dx ± idy). 
It is clear that, with the nondegenerate energy levels n' = (n,m), all of the previous arguments may be carried over 
directly to this example. 

In summary, we have shown a way to execute quantum control via Feshbach resonance - the density and also 
the quantum mechanical squeezing of a trapped atomic BEC may, in principle, be controlled by tuning the s-wave 
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scattering length. Experimental challenges such as the loss of atoms [21], and the inherent noise of the external 
magnetic field must be addressed in order for this scheme to be successful. 

This work was supported by NSF, ONR, and ARO. SC wishes to thank Prof. Keith Burnett for earlier discussions. 
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FIG. 2. Change in variance of the amplitude and phase quadratures [AXi(t)]^ and [AX2(t)]^ as a function of time. Time is 
given in units of 7r/a;t where u)t is the trap frequency. 
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